Shapefile

MN

# load data
eds_mn <- readOGR(dsn = "Data/eds", layer = "Manhattan_1910")

OGR data source with driver: ESRI Shapefile 
Source: "/Users/grace/Documents/QMSS/RA/19Fall/hnyc_street_dictionary/script/Data/eds", layer: "Manhattan_1910"
with 1480 features
It has 1 fields

lines_mn <- readOGR(dsn = "Data/lines", layer = "centerlines_MN_20190612")

OGR data source with driver: ESRI Shapefile 
Source: "/Users/grace/Documents/QMSS/RA/19Fall/hnyc_street_dictionary/script/Data/lines", layer: "centerlines_MN_20190612"
with 7919 features
It has 22 fields

# load map
tmap_mode("plot")

tmap mode set to plotting

tm_shape(eds_mn) + 
  tm_borders(col = "red", alpha = 1, lwd = 3, lty = 4) +
tm_shape(lines_mn) + 
  tm_lines(col = "blue", lwd = 1) 

BK

eds_bk <- readOGR(dsn = "Data/eds", layer = "Brooklyn_1910")

OGR data source with driver: ESRI Shapefile 
Source: "/Users/grace/Documents/QMSS/RA/19Fall/hnyc_street_dictionary/script/Data/eds", layer: "Brooklyn_1910"
with 1112 features
It has 1 fields

lines_bk <- readOGR(dsn = "Data/lines", layer = "centerlines_BK_20190612")

OGR data source with driver: ESRI Shapefile 
Source: "/Users/grace/Documents/QMSS/RA/19Fall/hnyc_street_dictionary/script/Data/lines", layer: "centerlines_BK_20190612"
with 1748 features
It has 20 fields

# load map
tmap_mode("plot")

tmap mode set to plotting

tm_shape(eds_bk) + 
  tm_borders(col = "red", alpha = 1, lwd = 3, lty = 4) +
tm_shape(lines_bk) + 
  tm_lines(col = "blue", lwd = 1)